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ABSTRACT 

We measure the Galactic rotation curve and its first two vertical derivatives in the first and fourth 
quadrants of the Milky Way using the 21 cm VGPS and SGPS. We find tangent velocities of the atomic 
gas as a function of galactic longitude and latitude by fitting an analytic line profile to the edges of 
the velocity profiles. The shape of the analytic profile depends only on the tangent velocity and the 
velocity dispersion of the gas. We use two complementary methods to analyze the tangent velocities: 
a global model to fit typical parameter values and a local fitting routine to examine spatial variations. 
We confirm the validity of our fitting routines by testing simple models. Both the global and local 
fits are consistent with a vertical falloff in the rotation curve of — 22±6 km s _1 kpc^ 1 within 100 pc 
of the Galactic midplane. The magnitude of the falloff is several times larger than what would be 
expected from the change in the potential alone, indicating some other physical process is important. 
The falloff we measure is consistent in magnitude with that measured in the halo gas of other galaxies. 

disk — Galaxy: fundamental parameters — Galaxy: 
— Galaxy: structure — radio lines: general 



Subject headings: ISM: general — Galaxy: 
kinematics and dynamics 



1. INTRODUCTION 

The shape of the Milky Way rotation curve is indicative 
not only of the integrated mass within a radius, but also 
of the local mass distribution and physical environment. 
The change in the rotation curve with distance from the 
plane can tell us about the gravitational potential of the 
disk as well as which physical processes affect the gas 
kinematics. Previous studies of H I near the plane as- 
sum ed that the rotation velocity was c onstant up to ~ 1 
kpc (|Celnik et al.lll979tlLockmanlll984 : in this paper we 
test that assumption. 

In other galaxies, the falloff of the halo gas rota- 
tion curve with distance from the plane is well es- 
tablished, having bee n observed in H I (NGC 891: 
Swaters et alll997lan df Oosterloo et al.ll2007l NGC 2403: 



Fraternaliet al.ll200l N GC 5775: iLee et al.ll200l Ha 



(N GC 891: iKamphuis et all [2007l NGC 5775: lRandll2000| 
andlTiillmann et al.ll2000f ) and X-rays (| Strickland et al.l 
|2004|). A recent study of the galactic halos of NGC 891, 
4302, and 5775 using H a found that the decrease in the 
rotation curve with scale height was a nearly constant 
— 15 to —25 km s" 1 per s cale height across the three 
galaxies (|Heald et al.ll2007h . In the Milky Way, obser- 
vations of an Ophiuchus superbubble in H I and H a at 
Galactic longitude ~ 30° and 3.4 kpc above the plane 
showed gas that lags 27 km s" 1 behind the rotation at 
6 = (jPidoprvhora et ail 12006). Halo gas motions are 
likely to be domi nated by complex physics, such as a 
Galactic fountain (|Shapiro fc FielJl976tlBregmanlll980f) 
where gas is pushed off of the disk, moves through the 
halo, and then returns to the disk. iBeniamir] (|2002l ) 
summarized several different mechanisms of establishing 
a vertical falloff in the rotation speed, including radial 
pressure gradients and magnetic tension. 

Nearer to the plane in the Milky Way H I disk, fluc- 
tuations in the velocity structure associated with spi- 
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ral arms were first observed by lOorfj ([1962). These 
were termed "rolling motions" because on one side of 
the plane there is a positive shift in the tangent ve- 
locity, and on the other side the sh i ft is t owards more 
negative velocities. lYuan fc Wallace! (|1973f ) tried to ex- 
plain the rolling motions as purely a geometric effect, 
but later work showed that this could not be true in 
all ca ses (|Strauss fe Poeppel Il976t iFeitzinger fc Spickerl 
119851 ). Rolling motions have also been observed in other 
spiral arm tracers, such as CO (jWouterlootJ Il98lh and 
stars (jBrosche fc Schwanl Il984f ). Note that it is diffi- 
cult to detect similar phenomena in the disks of external 
galaxies because ~ 10 pc spatial resolution and ~ 1 km 
s _1 velocity resolution are necessary. 

In principle, it is simple to extract the rotation struc- 
ture of the inner Galax y from the tangent velocities of gas 
at different longitudes (jKwee et al.l fl954) : for each longi- 
tude, just find the velocity where emission has fallen off. 
However, this task is complicated by several factors that 
obfuscate the determination of the tangent velocity, in- 
cluding velocity dispersion due to random motion of the 
gas, bulk flows often associated with spiral density wave s 
(|K^[l96l iShane fc Bieeer-Smithl fl96l lBurtonifl972l) . 
and noncircular motions, whether due to the Galactic 
bar (|Peterslll975r ) or some other cause. Nevertheless, H 
I spectra have been used to measure the inner Galaxy 
rotation curve for more than half a century. Early stud- 
ies also uncovered a large-scale asymmetry between in- 
ner Ga laxy rotatio n curves on either side of the Galactic 
center (lKerrl ll962fl. which may be a sign of a triaxial den- 
sity distribution iBlitz fc Spergel (|1991l ) . Much attention 
has been paid to bumps in the rotation curve as a func- 
tion of radius ( Shane fc Bieger-Sm ith 1966; B urtonlll971l : 
lAlvarez et allll990t iMcClure-Griffiths fc Dickevl l2007h . 
but comparatively little analysis has been done on the 
fluctuations above and below the plane. 

In this paper, we measure the inner Galaxy rotation 
structure close to the plane using 21 cm emission from H I 
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gas. For the first time, we allow the rotation curve to vary 
with height off the plane as well as with galactocentric 
radius. 

2. FORMALISM 

We generalize the existing rotation curve formalism to 
include variation in the rotation speed as a function of 
height off the plane. To accomplish this, we need two 
coordinate systems. The first is centered at the LSR, 
with coordinates Galactic longitude I, Galactic latitude b, 
and line-of-sight velocity v r . The second is a cylindrical 
coordinate system with the origin at the Galactic center, 
radial coordinate R, and height off the plane z. 

Analysis of the gas at the tangent points is straight- 
forward because, unlike the rest of the inner Galaxy, dis- 
tance determinations in these areas do not suffer from 
the distance ambiguity. Furthermore, simple geometric 
relations can be used to calculate the galactocentric ra- 
dius of a tangent point, R = Rq |sin£|, and the height 
of gas at a tangent point, z = i?o cos £ tan fr. Since this 
mapping is one-to-one at a tangent point, we use these 
two coordinate systems interchangeably. We assume Rq, 
the distance from the Sun to the Galactic center, is 8.5 
kpc. 

Assuming circular rotation, the relation for line-of- 
sight velocity is given by 



v r = sin(Z) cos(&) 



^-e(R,z)-e 



(i) 



where Oo is the rotation speed of the LSR (which we 
assume to be 220 km s _1 ), and we have generalized to 
a rotation curve that varies with height off the Galactic 
plane, <d(R,z). We assume cylindrical non-intersecting 
gas orbits, with no velocity in the z direction. At the 
tangent points, we then reduce eqn. ([I]) to an equation 
for the tangent velocity 



v t (l, b) = cos(6) [±9{R, z) - 9 sin(Z)] 



(2) 



where the "+" is for quadrant I (QI, 0° < I < 90°) and 
the "-" for quadrant IV (QIV, 270° < I < 360°). Wc 
use measurements of vt to constrain models of Q(R, z). 

3. DATA 

For the longitudes 18° < I < 67°, we used the 21cm 
H I V LA Galactic Plane Survey (VGPS) data (|Stil et all 
2006) to find the tangent velocities as a function of I 
and b. The VGPS is a combination of interferometric 
VLA observations with single-dish observations from the 
Green Bank telescope. The data have a velocity resolu- 
tion of 1.56 km s _1 , an angular resolution of 1', and an 
rms noise level of 2 K. The b coverage varies with longi- 
tude; near I ~ 20° the data are restricted to |6| < 1?3, 
while around I ~ 60° the range is |6| < 2?3. These data 
have a higher angular resolution than that of the Le i- 
den/ Argentine/Bonn survey, 0.5° (jKalberla et al.ll2005l) ; 
the better angular resolution provides the necessary res- 
olution in z. 

For the southern longitudes 270° < I < 358°, we 
used the 21 cm H I Southern Galactic P lane Survey 
(SGPS) data (|McClure-Griffiths et al.l[2005h . The SGPS 
is a combination of interferometric observations with the 
Australia Telescope Compact Array with single-dish ob- 
servations from the Parkes Radio telescope. These data 



have a velocity resolution of 0.8 km s _1 , an angular res- 
olution of ~2', and an rms noise level of ~1.6 K. The 
SGPS is limited to \b\ < 1?5, but to avoid edge effects 
associated with the interferometer coverage we further 
restricted to |b| < 1?4. 

To reduce these surveys to a more manageable size 
and minimize the impact of instrumental noise and small 
Galactic perturbations, we broke the data into 0?5 seg- 
ments in I and, for each segment, we constructed a com- 
posite spectrum by taking the median value of the bright- 
ness temperature at each velocity and for each b. This 
procedure reduced the number of data points in I while 
preserving the sampling in v r and b. We use this smaller 
data set for the remainder of this paper. 

Finally, we restrict to lines-of-sight with \z\ < 100 pc, 
to insure that each longitude has nearly uniform cover- 
age in z. We discuss this in more detail in £17.21 These 
cuts leave us with 37396 spectra in the north and 22897 
spectra in the south within the range 3 kpc < R < 8 kpc. 

4. TANGENT VELOCITY DETERMINATION 

4.1. Curve Fitting 

Various methods have been u sed to det e rmine v t for 
an H I spectra. For example, iMalhotral (|1995f ) chose 
the veloc i ty wh ere the profile crosses the 10 K thresh- 
old, iKerrl (|1962h chose the velocity of the first maximum 
in the spectra, IKerrl (|1964[ ) used the velocity where the 
brightness temperature cross es half the maximum, and 
iShane fc Bieger-Smithl (|1966l ) calculated v t using the in- 
tegral of the spectrum past the first maximum. Even the 
best of these methods requires correction, as modeling 
shows that the equivalent width velocity or the half max 
velocity differs from the actual t angent velocity by ~ 5 
km s" 1 dBurton fc Gordon! 1 1978h. A lso, functions such 
as Gaussians (IKulkarni fc Fichfl985h and error functions 
(|McClure-Griffiths fc Dickevl 1200711 have been fit to the 
emission falloff. In general, methods that rely on curve- 
fitting or the integral of the profile are more reliable than 
methods that depend solely on a temperature threshold 
or peak location because they will be less sensitive to in- 
strumental noise and perturbations in the H I gas density 
and velocity. 

Assuming that the disk ha s a uniform number density 
n and a spin temperature T dCelnik et al.1 (jl979t ) derived 
that the optical depth line profile close to the tangent 
velocity has the following form: 



t(v) = f{l,n,T s ,A) I exp 

10 



1 



[w 



dw (3) 



where t = (v — v t )/a, a is the velocity dispersion, / 
is a function that parameterizes the amplitude of the 
line, and A is the Oort constant evaluated at the tan- 
gent point. The shape of the profile depends only on the 
velocity dispersion and the tangent velocity of the gas; 
other parameters such as the density and shape of the 
rotation curve affect the amplitude. We convert from r 
to brightness temperature Tb by assuming the spin tem- 
perature is 155 K everywhere, larger than the hi ghest Tb 
in our data set. iRohlfs fc Kreitschmannl (|1987f ) further 
refined eqn. by allowing the kinematics of the gas to 
deviate from circular symmetry and adding in radial mo- 
tions; for the purposes of this paper, such a high level of 
detail is not necessary. 



Milky Way Rotation Curve and Vertical Derivatives 



3 



The irregular shape of many profiles sometimes makes 
it difficul t to choose which peak sho uld be fit with 
eqn. ©; Roh lfs fc Kreitschmannl (|1987l ) were forced to 
fit spectra interactively, determining the limits for each 
fit by eye. We use a more objective technique: first, we 
smooth each spectra with a 4 km s _1 boxcar window. 
Next, we find the location of the first local maximum in 
the smoothed spectra with a brightness temp larger than 
10 K; we call this velocity v max . We then fit the origi- 
nal spectrum with eqn. ([3]) using a Levenberg-Marquardt 
non-linear least squares fitting routine with three free pa- 
rameters: v t , amplitude, and a. We include in this fit 
data from velocities in the range w max — 8 km s _1 < v < 
t) mm + 65 km s _1 for QI and v max — 65 km s _1 < v < 
^max + 8 km s _1 for Q1V. These limits were optimized 
by minimizing the errors in parameters determined from 
our models, described below in $721 

The fits are generally a good match to the emission 
spectra. However, there are some exceptions. Some spec- 
tra are too irregular for the fitting routine to converge on 
a physically reasonable solution; an intermediate velocity 
cloud can fall in the range of fitted velocities, or there can 
be a particularly large bulk flow that distorts the shape 
of the spectrum. In cases like these, v t found by the 
fit often falls outside the velocity boundaries of the fit. 
When this happens, we flag the spectra at that (I, b) and 
remove its Vt from the rest of the analysis. 0.6% of the 
spectra in QI and 0.7% of the spectra in QIV are flagged 
in this way. 

The map of fitted Vt(R, z) is shown in Fig.[T]for QI and 
Fig. [2] for QIV. Maps on both sides of the Sun-Galactic 
center line show the expected large scale trend of lower 
\vt\ at larger R. The Vt surfaces are fairly smooth with 
a few small pockmarks. In a system with no change in 
the rotation curve with z, the contours in these maps 
would be nearly parallel to the z axis, since cos (6) ~ 1. 
This is clearly not the case in many regions; for one ex- 
ample, examine the —90 km s _1 contour in the southern 
map, which shows evidence for a dvt/dz term. Many 
contours in the \v t \ = 50 km s _1 to 90 km s _1 range 
appear several hundred pc closer to the Galactic center 
in QIV compared to QI (see Fig. [3|) . This could result 
from a genu ine asymmetry in the QI/QIV rotation curves 
(|KerHfl962l ). deviations from the axial symmetry of the 
potential, or from bulk flows and other peculiar motions. 

In Figs. 0] and [5l we plot the tangent velocity on top of 
Tf, contours for several I. Note that the tangent velocity 
points do not simply follow a single brightness tempera- 
ture contour; in general, they are at lower Tf, for larger 
z, because the overall amplitude of the spectra decreases 
away from the plane as the density falls off. On the ve- 
locity scale of these figures the Vt are quite regular, albeit 
with a few outliers. 



4.2. Removing v t Surface Discontinuities 

The tangent velocity surface should be a smooth func- 
tion of I and b, but the surfaces derived in S I4.1I have 
discontinuous bumps. For example, local maxima in the 
spectra can result in large discontinuities in v t {z), as in 
the / = 335° panel of Fig. [5] near z ~ 70 pc. We flagged 
discontinuities like this one using a median-type filter and 
did not include them in our analysis. The filter worked 
as follows. 



1. Calculate the function 

M(k, bj) = Median [v t (l l+k , b j+m )] 



(4) 



where — 1 < k < 1 and —60 < m < 60 for QI and 
—30 < m < 30 for QIV. In other words, M is the median 
of all the points inside a rectangular box in (7, b) space 
with a width of 3x121 elements for QI and 3x61 elements 
for QIV; the allowed ranges of k and m set the size of the 
median filter in the / and b dimensions, respectively. The 
box encompassed more elements in b than in I because we 
already reduced the number of elements in I as discussed 
in Sj3] Twice as many points were included in b in QI 
because of the higher resolution of the VGPS compared 
to the SGPS. The stepsize in / for both quadrants was 
0?5, while in b it was 0?005 for QI and 0?011 for QIV. 
Points that were flagged due to poor fits in E14.ll were not 
included in the median calculation. 

2. If \v t (li,b j )-M(l i ,b j )\ > 5km s" 1 then flag v t (h,bj) 
and do not include it in our analysis. This threshold 
eliminated most of the bumps that were inconsistent in 
adjacent longitudes. 

Flagged data points for a few I are shown with white 
x's in QI (Fig. [4]) and QIV (Fig. [5]). 4.2% of the tangent 
velocities in QI were flagged by this filter, as were 6.9% 
of the velocities in QIV. This filter generally flagged the 
Vt points where the tangent velocity varies erratically. 

Discontinuities that persist over a large angular scale 
were not removed by this filter. In particular, a large 
discontinuity is located between I — 57° and 60?5 with 
-50 pc < z < -100 pc (see Fig.Ej I = 60°). We flagged 
this region by hand and removed it from the analysis. 

5. GLOBAL ROTATION CURVE FIT 

To model the global rotation curve, we extende d a 
linear rotation curve explored in iFich et ail (|1989h by 
adding two vertical derivative terms: 



0(R,z) = ai e Q + a 2 ^R- 
Ro 



d& . d<d 

M^ + lTz z - 



(5) 



The global rolling motion is given by dQ/dz, and the 
global falloff from the plane by dO/d\z\. Throughout 
this paper, we call these the rolling parameter and the 
falloff parameter, respectively. Using eqn. ([2]), we then 
have an equation for the observed tangent velocities: 



Vt(l, b) — cos(6) 



±ai8 ± a 2 ^-R 
Ho 

d& dQ 

± TTTT \ z \ ± ~T Z ~ @osm(/) 
d\z\ dz 



(6) 



where the "+" is for QI and the "-" for QIV. Alter- 
natively, we could have replaced the third term with 
z 2 d 2 Q/dz 2 , but this choice did not affect our results. 

Since points are unequally spaced in (i?, z) space, per- 
forming a simple linear least squares fit would unevenly 
weight some R more than others. Also, the stepsize in b 
in QI is smaller than in QIV, so QI would have a larger 
influence on a combined fit. Furthermore, although there 
are more than a hundred spectra at each longitude, many 
of the derived tangent velocities are not independent. 
Bulk flows and other peculiar motions are responsible 
for most of the deviations from circular motion; these 
motions have a typical velocity scale of 8 km s . Small 
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Fig. 1. — Tangent velocity Vt as a function of R and z in QI. Color contours are drawn every 10 km s — 1 and are marked on the color bar 
on the right in km s — 1 . The data beyond R = 7.8 kpc are missing due to the / limitations of the VGPS. Dotted lines encircle regions that 
have been flagged, either by the vt fitting routine or by the smoothness filter; these regions should be ignored. The —40 km s^ 1 and —90 
km s — 1 contours are marked with solid black lines. 
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Fig. 2 
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— Same as Fig. [T] but for QIV. The missing data with R > 7.5 kpc and \b\ > 1?4 is due to the restriction on b for the SGPS. The 
1 and 90 km s _1 contours are marked with solid black lines. 



regions can vary from the local rotational speed by 20- 
30 km s" 1 (|Brand fc Blitz! Il993f ). The physical scale of 
these bulk flows is difficult to determine, but we estimate 
there are approximately f independent data points with 
\z\ < 100 pc for each /. 

We construct an interpolated data set with a z grid 
given by —90, —70, —50, . . . , 90 pc but leaving unchanged 
the spacing in I. We flagged points in the interpolated 
grid if cither of the two points adjacent in b in the original 
grid were flagged. We performed a linear least squares fit 
to the interpolated data within the range 3 kpc < R < 8 
kpc and \z\ < 100 pc. The results for the Ql, QIV, and 
combined data set are shown in Table [T] along with the 
formal errors in the fit. We emphasize that these rotation 
curve parameters are valid only inside \z\ < 100 pc; we 
have not constrained the falloff outside of this region. 

What is the effect of the bulk flows on the values of 
the fitted parameters? In short, quantifying the impact 



of bulk flows requires a self-consistent simulation of gas 
dynamics in the Milky Way disk, which is beyond the 
scope of this investigation. We have assumed that these 
random motions do not have a strong effect on the values 
of the fitted parameters, other than to correlate nearby 
data points. 

Are there places in the inner Galaxy where this anal- 
ysis is inappropriate? Noncircular orbits can result from 
the presence of a stellar bar; this would invalidate an 
assumption that led to eqn. ([2]). A recent survey of 
mid-infrared sources toward the inner Galaxy detected 
a QI-QIV asymmetry in star counts fo r I < 30°; this 
corre sponds to a bar radius of ~4.4 kpc (|Beniamin et alj 
l2005h . We also performed a global fit restricted to spec- 
tra with 4.4 kpc < R < 8 kpc, but this did not strongly 
affect the derived parameters. We defer an in depth dis- 
cussion of the physical interpretation of the parameter 
values until 36l 
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Fig. 3.— Contour lines for \v t \ = 30,50,70,90,110 km s" 1 for QI (solid) and QIV (dotted). The contour values are given by the color 
bar. 



TABLE 1 

Global Least-Squares Fit Parameters 



Data 


ai 


U2 


de/d\z\ 


dO/d\z\ COII 


dS/dz 


QI 
QIV 
Combined 


0.855±0.004 
0.829±0.004 
0.844±0.003 


0.209±0.006 
0.221±0.005 
0.212±0.004 


-19±7 
-12±7 
-16±5 


-26±8 
-16±8 
-22±6 


-8±4 
9±3 
1±2 



Note. — The global rotation curve assumes Ro = 8.5 kpc and 0o = 220 
km s _1 . The parameters a\ and ci2 are unitless, d&/d\z\ and dQ/dz are 
in km s — 1 kpc -1 , and d& /d\z \ CO rr is the falloff parameter corrected for the 
systematic error discussed in q7.2. II 



6. LOCAL ROTATION CURVE 
PARAMETER DETERMINATIONS 

Do the vertical derivatives of the rotation curve vary as 
a function of Rl We investigated by fitting the vertical 
derivatives at each I independently. We approximated 
the rotation curve by the first three terms of its Taylor 
expansion, 



e(R, z )*e(R,Q) + \z\^ 

d\z\ 



de 

dz 



(7) 



We then wrote the equation for the tangent velocity in 
terms of the 6 = solution to eqn. @, Vtfl. 



v t (b) = cos(6) 



Vtfl ± 



de 



de 




±z ^r 




2=0 dz 


z=0 



(8) 



where again the "+" is for QI and the "— " for QIV. 
Note that this equation has the observables Vt(b) and 
three unknowns: v t .o, d@/d\z\, and de/dz. For each I 
in our sample, we determined the three rotation curve 
parameters in eqn. ([5]) using a linear least squares fit. 
Several fit samples in QI (Fig. QI and QIV (Fig. [5]) are 
shown as solid lines on top of the contours and v t fits. 

Again we have the problem of bulk flows correlating 
the tangent velocities at nearby points. Ideally, we would 
proceed by calculating the correlation matrix to use in 
the fit, but the physical scale of the flows is not well deter- 
mined. Instead, we adopt an ad hoc solution: calculate 
the formal error in the fit assuming the data points are 
independent, and then scale up the errors by assuming 



there are only ten independent data points between —100 
pc < z < 100 pc, i.e. we multiplied the formal errors by 
\J number of points/10. 
We now discuss each of the three fitted parameters. 

6.1. Rotation Curve at z = 

In Figs. [5] and [7] we plot the parameter 0(i?, 0) in QI 
and QIV, respectively. Note that, for the purposes of 
this plot, we have assumed Oo = 220 km s _1 with no 
error, but the determination of v t does not depend on 
the value of 6 . The SGPS for R > 7.5 kpc does not 
cover the full range of \z\ < 100 pc because of the lim- 
ited survey extent in b. For these longitudes, we remove 
the fits of the two derivatives from our analysis, but keep 
the easier to determine 0(i?, 0). In both QI and QIV the 
rotation curve rises as R increases, consistent with pre- 
vious measuremen ts of the inner Galaxy rotation curve 
(|Kwee et al.lll954f ). On top of the O(i?,0) fit points, we 
plot the global fit to the rotation curve from Sj5] The fit 
traces the measured rotation curve well, reinforcing the 
conclusions of our global analysis. 

Our fit to the QIV rotation curve is approxi- 
mately 7 km s -1 lower than th e one determined in 
iMcClure-Griffiths fc Dickevl (|2007j ). We suspect this dis- 
crepancy is due to the differences in vt fitting methods. 
In both q uadrants, our data a re not well matched by the 
fit from iBrand fc Blit3 (|1993l ). However, data fr om the 
outer Galaxy was also used to constrain their fit, so this 
is not surprising. 
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25.0 101.9 46.0 31.8 30.0 108.7 -113.7 -22.3 




-100 -50 50 100 -100 -50 50 100 
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Fig. 4. — The measured brightness temperature as a function of z and v r for several QI I. Contours are drawn every 10 K and lighter 
contours denote higher Tj,. The x's mark vt for each b, and the fit to the tangent velocities is shown with a solid line. Flagged tangent 
velocities are marked with white x's. The fitted parameters for these I from [0are written on top of each plot in the following order: I 
(degrees), Vt (km s —1 ), dO/d\z\ (km s -1 kpc - 1 ), and dS/dz (km s -1 kpc - 1 ). 
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Galactic Longitude (deg) 




Fig. 6. — The rotation curve Q(R, 0) as measured at the tan- 
gent points in QI. The solid line is the rotation curve fit to the 
QI data from £|5l an d the dashed line is the fit to QI data from 
IBrand fc Blit3(|l993h . 



Gclactic Longitude (deg) 
325 320 




R (kpc) 

Fig. 7.— Same as Fig. [6] but for QIV. The solid line is the 
rotation curve fit to the QIV data f r om j|5l the dashed line is the 
fit to QIV data from lBrand fc Blitd 111993). and the dotted line is 
a fit to the SGPS data from lMcClure-Griffiths fc Dickevl IpfflOW ). 

The QI and QIV rotation curves are shown on the same 
plot in Fig. [51 The asymmetry in the northern and south- 
ern rotat ion curves i s quite similar in shape to that mea- 
sured bv lKerrl (|1964h . who observed that the QI curve is 
larger than the QIV curve between 4.2 kpc < R < 6.8 
kpc. Compared to this previous work, the magnitude of 
our curves is different due to the change in the preferred 
values for Galactic co nstants. More recent work shows 
the same asymmetry (jBlitz fc Spergellfl991l ). 

Bumps in the inner Galaxy rotation curve have long 
been interpreted a s the locations where spir al arms cross 
the tangent circle (|Yuanlll969l: lBurtodll97l|) . We will not 
speculate on the meaning of the bumps, other than to say 
they are about 10 km s -1 in amplitude , as expected from 
density wave theory (see lBurtonlll97ll ). 

6.2. Falloff Parameter 




R (kpc) 



Fig. 8. — The rotation curve @(R, 0) as measured at the tangent 
points in the inner Galaxy. The x's connected by a solid line are 
the QI values; the triangles connected with a dashed line are the 
QIV ones. The rotation curve from the fit to the QI and QIV data 
(Table [TJ is shown with the heavy solid line, the QI fit is shown by 
the light grey line, and the QIV fit is shown by the dark grey line. 

The falloff parameter dQ/d\z\ is shown for QI in Fig. [9] 
and QIV in Fig. [lOl In QI, the falloff appears centered 
near —30 km s -1 kpc -1 . We include a correction for a 
~ 5 km s _1 kpc -1 systematic error due to the falloff of 
gas density with z discussed in §7.2.21 The falloff pa- 
rameter is roughly constant with radius, although small 
oscillations are apparent throughout, and a large jump 
is seen near R ~ 4.5 kpc. In comparison, the falloff in 
QIV is much less regular. The points are also centered 
around —30 km s -1 kpc" 1 , but the scatter is much larger 
and their values are less correlated from one point to the 
next. There also is a large group of points near R ~ 7 
kpc where d@/d\z\ > km s -1 kpc -1 ; it is unclear what 
the significance of this group of points is, and also why 
the QI and QIV plots have such a different appearance 
in detail. 

Overall, the local analysis again confirms the global 
fit because they both point towards a falloff parameter 
with magnitude approximately —20 km s _1 kpc -1 . The 
roughly constant value of d<d/d\z\ in QI consistent with 
the global fit from pleads us to conclude that the falloff 
is a global feature of the inner Galaxy rotation curve (at 
least within 100 pc of the plane), and is not limited to 
any one region. 

6.3. Rolling Parameter 

The rolling parameter dO j dz is shown for QI in Fig. 1111 
and QIV in Fig. [T2] In both these figures, nearby points 
are highly correlated and the curves have several peaks 
and valleys. Both plots appear to be centered near zero 
velocity. 

There is an asymmetry between the QI and QIV val- 
ues of the rolling parameter in that the QIV data have 
higher amplitude oscillations than those in QI. The high 
level of rolling motion in the QIV tangent circle may be 
the reason the QIV falloff parameter has high scatter 
compared to that of QI. Though the two derivatives are 
nearly orthogonal, the falloff may be easier to detect in 
a region with quiescent H I and the QIV tangent circle 
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Fig. 9. — The radial dependence of d&/d\z\ as measured at the Fig. 11. — The radial dependence of d@/dz as measured at the 

tangent points in QI. Error bars are from the adjusted error in the tangent points in QI. Error bars are from the adjusted error in the 



fit (36}. 



fit (<JB). 
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Fig. 10.— Same as Fig.|9]but for QIV. 



Fig. 12.— Same as Fig.[TTJbut for QIV. 



clearly has a good deal of kinetic energy due to rolling 
motions. The correlation length of the data in QIV also 
appears to be larger than that of QI. 

It is difficult to compare the rollin g motions in these fig- 
ures to those observed previously bv lFeitzinger fc Spickerl 
(1985). The authors claim typical values of ±20 km s _1 
kpc -1 and values reaching —75 km s _1 kpc -1 in localized 
regions. They determined these numbers by tracing the 
tilt in Tfc contours specifically in spiral arms, which are 
not directly comparable to our calculations. Neverthe- 
less, the two methods of measuring rolling motions sat- 
isfyingly produce slopes of the same order-of-magnitude. 

Because of the oscillatory nature of the rolling motions 
in both quadrants, as well as the global fit consistent with 
zero slope, wc conclude that the rolling motions average 
to zero over the inner Galaxy (again, within 100 pc of the 
plane). However, local fits demonstrate that slopes as 
large as ±100 km s -1 kpc -1 appear in localized regions. 



7. FIT VARIATIONS AND SYSTEMATIC ERRORS 
7.1. Fit Variations 

How dependent are the parameter values on the 
method of fitting v {l We refit all the da t a usin g the 
integral method of I Shane &; Bieger- Smith! ( |1966l ) , and 
find that the vt are generally ~5-10 km s _1 larger in 
magnitude than those d etermined with the curve fitting 
method, as predicted bv lBurton fc Gordon! (|1978f ). This 
results in a rotation curve with higher velocities. If wc 
run this data set through our global fitting routine, the 
falloff parameter for the QI and QIV combined data is 
— 12 ± 4 k m s" 1 kpc -1 , withi n la of the result from Sj5] 
using the ICelnik et al.1 (|1979l ) vt fitting method. The 
rolling parameter has a global value of 1 ± 2 km s _1 
kpc -1 , essentially a match. We conclude that our use 
of the analytic model fitting method does not bias the 
appearance of the vertical derivatives. 

We have also tried replacing z in eqn. {8j with (z— z m ), 
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where z m is the center of a Gaussian fit to the brightness 
temperature as a function of z for a given R (in other 
words, the midplane of the gas distribution). The mid- 
plane can be displaced by up to 50 pc from z — in 
the inner Galax y, though smaller values are more typical 
(|Malhotralll995ft . We ran the local fit to the QI data with 
this modification, and saw little change in the values of 
the two derivatives. 

The Sun is located about 1 5 pc a bove the Galac- 
tic pla ne (|Freudenreich et ail 11994] and references 
therein. ICohenlll995l lNg et al.lll997l iBinnev et~aTlll997l 
iMagnani et alJ Tl985h . We checked if this influenced 
our results by recalculating the vertical derivatives 
with an adjusted z coordinate, z a dj = z + 15 pc — 
i?o cos(7) tan(0.1°). The change in ai, a2, and d<d/d\z\ 
were all within the la errors, but the combined value for 
dQ/dz became 6 ± 2 km -1 kpc -1 . 

7.2. Systematic Errors 
7.2.1. Testing Global Fit Accuracy with Models 

For the purpose of testing our global rotation curve 
parameter determinations (©, we created simple models 
of Galactic H I emission with the same I and b spacing 
and v r resolution as the SGPS. We assumed a constant H 
I surface density of 4 M Q pc -2 (|Burton fc Gordon|[l978ft 
and a two-component gas, 80% with a velocity dispersion 
of 8 km s -1 and the remainder with a dispersion of 35 
km s" 1 (|Kulkarni fc Fichl [l985h . The gas in the model 
is ass umed to follow a sech 2 distr ibution in z (ISpitzerl 
11942ft with a scale height of 120 pc (|Lockmar]ll984ft . and 
for simplicity we did not change the proportions of the 
different velocity dispersion components as a function of 
height. We included optical depth effects using a spin 
temperature of 155 K. We constructed several different 
models for the change in the rotation curve as a function 
of R and z (Table [2|) . We did not add random noise 
to the models because we are interested in estimating 
the systematic uncertainty solely due to the falloff in gas 
density with \z\. 

After reducing the size of the model data using the me- 
dian method described in <J3]and fitting the v t as in 
we performed the linear least squares fit described in §3] 
to recover the parameters of Q(R, z) in each case. Table 
[2] shows that the input parameters were essentially recov- 
ered, suggesting a small systematic error in all parame- 
ters except dQ/d\z\, where the error is generally 4 ± 1 
km s -1 kpc -1 . This error occurs because we fit v t at 
latitudes where b ^ 0, and the assumption in eqn. 
that n is constant along the line of sight is not justified. 
On average, points further along the line of sight will 
have smaller n than more nearby points, since the far 
points are at a higher z and the gas is generally densest 
at b = 0. This varying gas density causes small inac- 
curacies in the fitted v t which translate to a systematic 
error in the falloff parameter. The surface density of 
H I plays a role in setting the magnitude of this effect. 
We correct for this systematic error in the column la- 
beled d9/d|z| corr in Table [T] The systematic errors in 
a\ and ai are sometimes significant, but we do not cor- 
rect for them here. We also modeled QI with a\ = 0.82, 
a 2 = 0.22, dQ/d\z\ = -30 km s -1 kpc -1 , and dQ/dz = 0, 
similar to the measurements in $5j and we found the er- 
ror in the falloff parameter was 7 km s _1 kpc -1 . When 



the QI and QIV data sets were combined, the error in 
the falloff parameter was 6 km s -1 kpc -1 . From our ex- 
perience with the QIV models, we assigned an error in 
these systematic corrections of ±1 km s _1 kpc -1 . 

7.2.2. Testing Local Fit Accuracy with Models 

To test our local parameter determinations (Sj6j), we 
created simple models of Galactic H I emission at I — 
317?5 (i?=6.2 kpc) with the same I and b spacing and 
v r resolution as the SGPS. We used the same gas char- 
acteristics as our previous models in §7.2.11 except we 
tried several different velocity dispersion combinations. 
The basic model has a two component gas, 80% with a 
velocity dispersion of 8 km s -1 and the remainder with 
a dispersion of 35 km s" 1 (jKulkarni fc Fichl fl985h . A 
recent analysis has suggested a three component model 
with 56% of the gas with a velocity dispersion of 6 km 
s" 1 , 38% with 12 km s" 1 , and 6% with 26 km s -1 
gas (jMcClure-Griffiths fc Dickevl 120071 ): we construct a 
model with these parameters as well. For completeness, 
we also constr uct a model made solely of 8 km s -1 gas 
(|Burtonlll976ft . We made models with various values of 
velocity dispersions and vertical derivatives. 

We performed the fit described in §4.11 to determine 
v t (l = 317?5,6) for these models, and tried to recover 
the parameters of Q(R = 6.2 kpc, z) in each case (Table 
[3]). The input parameters are essentially recovered, sug- 
gesting systematic errors of ~ 1 km s -1 in vt,o and ~ 5 
km s -1 kpc -1 in the falloff again due to the change in gas 
density with z. For models with more than one velocity 
dispersion, it is possible to perform a more complicated 
fit by summing multiple versions of eqn. ([3]), one for each 
component, but we found that this does not improve the 
accuracy of the fitted values of O and its derivatives. 

Although the systematic errors in the falloff parameter 
local analysis do not vary as a function of the gas velocity 
dispersion, they do vary as a function of R. To correct 
for this effect, we used the global model with parameters 
ai = 0.82, a 2 = 0.22, d@/d\z\ = -30 km s -1 kpc -1 , and 
dQ/dz = km s -1 kpc -1 and ran it through our local 
analysis. The routine returned errors in dQ/d\z\ ranging 
from ~ —1 km s _1 kpc -1 at R = 3 kpc to ~ —12 km s -1 
kpc -1 at R — 8 kpc. These magnitude of these errors 
is expected from the errors in the global analysis. We 
corrected the local measurements of the falloff parameter 
in Figs, land [101 

8. DISCUSSION 

Averaging spectra over b, as other authors have done 
(|Malhotral 119951 : iMcClure-Griffiths fc Dickevl 12007ft . can 
result in errors of ~l-2 km s -1 in the determination of 
Q(R, 0) because the rotation speed changes with height. 
This effect is more harmful in surveys that are limited 
in b (compared to a survey limited in z), because their 
range of z varies as a function of I. A systematic error 
can therefore be introduced. A similar problem occurs if 
the data do not extend to the same distance in b on either 
side of z — 0. If it is necessary that spectra be averaged 
to improve signal to noise, we recommend averaging over 
a set limited in z rather than one limited in b. 

Is the magnitude of the falloff parameter surprising? 
What change in the rotation speed is expected due the 
shape of the gravitational potential? We assumed the po- 
tential is dominated by the stellar disk contribution. At 
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TABLE 2 
Global Model and Fit Parameters 



01 


02 


d@/d\z 


Model 
dQ/dz 


CT1 


C2 0-3 


Aai 


Aa 2 


Fit 

AdO/d\z\ 


Ad&/dz 


1.00 


0.00 


-30 


-30 


8 


35 


-0.009 


0.007 


4 


1 


1.00 


0.00 


30 


30 


8 


35 


-0.009 


0.007 


4 


1 


0.82 


0.22 


-30 


30 


8 


35 


-0.014 


0.013 


5 


1 


0.82 


0.22 


-30 


-30 


8 


35 


-0.014 


0.013 


5 


1 


0.82 


0.22 


30 





8 


35 


-0.013 


0.012 


4 





0.82 


0.22 








8 


35 


-0.013 


0.012 


4 





0.82 


0.22 


-30 





8 


35 


-0.014 


0.013 


5 






Note. — The derivatives d&/d\z\ and dQ/dz are given in km s kpc . The param- 
eters ai and ai from eqn. are unitless. Velocity dispersions a are given in km s . 
The systematic errors given in the fit column are defined as Ai = XR t — x mo< j e l- 



TABLE 3 
Local Model and Fit Parameters 



Vt,0 


dB/d\z\ 


Model 
dQ/dz 


°"i 


o~2 


0-3 


Av t ,o 


Fit 

AdQ/d\z\ 


AdQ/dz 


-71 


30 





8 


35 




1 


5 





-71 








8 


35 




1 


5 





-71 


-30 





8 


35 




1 


5 





-71 


-60 





8 


35 




1 


5 





-71 


-90 





8 


35 




1 


5 





-71 


-30 





8 






1 








-71 


-30 





6 


12 


26 


1 


2 





-71 


-30 


-30 


8 


35 




1 


5 





-71 


-30 


30 


8 


35 




1 


5 






Note. — Tangent velocities vtfi and velocity dispersions a are given in 
km s — 1 . The derivatives dQ/d\z\ and dQ/dz are given in km s — 1 kpc -1 . The 
systematic errors given in the fit column are defined as Ax = Xfl t — x mo< j e i- 



iPidoprvhora et all (|2006f ) measured a falloff of 27 km s _1 
for a single structure 3.4 kpc above the plane; the falloff 
derived from this object is not directly comparable to our 
measurement. Since we measure the rolling motions and 
falloff parameter with \z\ < 100 pc where the derivatives 
have magnitudes of < 10 km s -1 , our result will also 
be difficult to confirm in o ther galaxies. At higher \z\, 
iFraternali fc Binneyl (|2006f ) examine the halo gas kine- 
matics in NGC 891 and NGC 2403 and conclude that 
the loss of angular momentum is due to the interaction 
between accreted low-angular momentum extragalactic 
gas and halo gas. It is conceivable that this effect is tied 
to the falloff we observe in the atomic disk gas, possibly 
with the vertical magnetic field playing a role. Addition- 
ally, it may be important to consider rolling motions and 
the falloff parameter when examining gas kinematic phe- 
nomena such as supernova remnants and H II regions; 
the rotation curve derivatives can shear shells of gas if 
they are large enough. 

9. SUMMARY AND CONCLUSIONS 

We fit an analytic form to the VGPS and SGPS H 
I spectra and found the tangent velocity as a function 
of Galactic longitude and latitude. Using these tangent 
velocities, we fit a simple global model for the rotation 
curve and its first two vertical derivatives for lines of 
sight 100 pc or less from the plane. We then determined 
the rotation speed, rolling motion, and falloff in rotation 
speed from the plane as a function of Galactic radius. We 
tested the accuracy of our fitting routines using a variety 
of simple H I model spectra. Our local parameter values 
are in agreement with our global fits. We demonstrated 



some locations in the disk this assumption will not hold, 
but spherical components like the bulge and the halo re- 
sult in negligibly small dQ/d\z\ at \z\ = 100 pc. Next, we 
assumed that the gas travels in non-intersecting cylindri- 
cal orbits, with the rotation speed exactly balancing the 
radially inward gravitational force, gn, resulting in the 
relation Q(R,z) — \/—RgR . For an exponen tial vertical 
distribution of stellar mass. lPeng et al.1 (|2002f ) found that 

l-OO H e ~ a \ z \ _ g— 

g R (R,z) = -2irG / pS Q {p)J 1 {R0) a d(3 

J ° (I) -1 

(9) 

where 

So(fi) = / rJ Q ((5r)^{r)dr. (10) 
Jo 

Jo and Ji are the zeroth and first order Bessel functions, 
is the stellar surface density, and a is the inverse 
of the stellar scale height h. We solved these equations 
using h = 335 pc, E(R) = S exp(-R/R d ), R d = 2.5 
kpc, and S(i? ) = 54 M Q pc~ 2 . Within the range 3 
kpc < R < 8 kpc, this model predicts |ef0/d|z|| < 5 km 
s _1 kpc -1 , too small to explain the magnitude of the 
observed effect by about a factor of four. Some other 
process besides simple gravitational physics must be at 
work. 

In Table 0] we summarize measurements of the falloff 
in other galaxies. Though these measurements rely on 
emission from halo gas kiloparsecs from the plane, the 
falloff parameter we measure is still remarkably con- 
sistent in magnitude. Additionally, in the Milky Way 
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TABLE 4 

Falloff Parameter Measurements in Galactic 
Halos 



Galaxy 


Tracer 


d&/d\z\ 




Reference 


NGC 891 


H I 


-25 to - 


100 


Swaters et al. (1997) 


NGC 891 


H a 


-15 to - 


-18 


Heald et al. (2006a) 


NGC 891 


H I 


-15 




Ooaterloo et al. (20071 


NGC 891 


H a 


-18.8 




Kamphuis et al. MiST) 


NGC 2403 


H I 


-25 to - 


-50 


Fraternali et al. (20021 


NGC 4302 


H a, N II, S II 


-30 




Heald et al. (20077" 


NGC 5775 


H a 


-8 




Heald et al. (2006b) 


Note. - 


The derivative dQ/d\z\ is g 


iven 


in km s" 1 



kpc 



that the rolling motions are local effects, while the falloff 
from the plane has a consistent value of —22 ± 6 km s _1 
kpc -1 , too large to be explained by gravitational physics 
alone, but consistent with the falloff measured in the 
halos of other galaxies. 
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